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New Matter Effects and BBN Constraints for Mass Varying Neutrinos 
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The presence of light (tua ~ 10 _6 eV) scalar fields in the early universe can modify the cosmology 
of neutrinos considerably by allowing their masses to vary on cosmological times. In this paper, we 
consider the effect of Planck-suppressed couplings of this scalar to electrons and show that such cou- 
plings can easily make new sterile states thermally inaccessible in the early universe, preserving the 
successes of big bang nucleosynthesis predictions. We consider the circumstances under which these 
effects give the proper initial conditions for recently considered models of neutrino dark energy, and 
consider limits from tests of the equivalence principle. The parameters which satisfy cosmological 
constraints naturally give rise to interesting signals in terrestrial neutrino oscillation experiments. 



I. INTRODUCTION 

One of the most exciting discoveries in recent years has 
been the demonstration of the existence of neutrino os- 
cillations. The angular flavor dependence of atmospheric 
//-neutrinos pj gave solid evidence to the phenomenon, 
which was followed up by controlled terrestrial experi- 
ments, such as K2K 0] and KamLAND Q. This, to- 
gether with the SNO results 0, has now solidly estab- 
lished neutrino oscillations as the explanation of the long 
standing solar neutrino problem[5|, la LZl la l2| • 

However, many questions remain. Namely, what is the 
origin of neutrino mass and what sets the scale? Are 
they Dirac or are they Majorana in nature? What is the 
impact of neutrino mass on cosmology? 

One exciting possibility is that neutrino mass is associ- 
ated with a new, light, scalar field. This is quite reason- 
able, given that other fundamental fermions have mass 
arising from the Higgs boson, and that radiative correc- 
tions to the mass of such a scalar can be controlled [To| . 
New scalar forces for neutrinos have been considered for 
many years [Til Il2| . It has been shown that such a field 
could be the source of all neutrino mass within matter 
[l3|. Additionally, others have considered the possibil- 
ity that a slowly rolling scalar field with mass 10 _33 eV 
might change the mass of neutrinos on Hubble times |l4j | , 
possibly impacting questions of leptogenesis 0] . 

Recently, it has been proposed that the dark energy 
might arise due to interactions between relic neutri- 
nos 0| ■ Here, no Hubble mass particles need be in- 
voked, but rather particles with masses in the range 
10 -6 eV— 10 -8 eV. More recently, supersymmetric mod- 
els have been presented as a means to provide theories 
where all masses and couplings take on natural values 
O- (See also [13 for an alternative supersymmetric 
model.) One exciting possibility is that these new forces 
might allow us to test dark energy in neutrino oscillation 
experiments arising from medium dependent mass val- 
ues [3 • Studies show that there is still a great deal of 
room within neutrino experiments for such effects |19| . 
with some studies argu ing for possibly interesting effects 
on solar neutrinos yQ, |2J . 



A. The early universe problem of sterile neutrinos 

Models of this sort almost always involve the presence 
of a singlet neutrino, which, at present times, for natu- 
ralness reasons, has a mass of O(eV). The introduction 
of such a particle forces us to consider the early universe 
behavior of the theory, as such a neutrino is generally 
populated in the early universe, leading to severe con- 
straints from big bang nucleosynthesis [22|. 

Moreover, the neutrino models of dark energy require 
that the relic neutrinos be mass eigenstates, not interac- 
tion eigenstates, which raises another question of initial 
conditions. 

In this letter, we shall see that the the presence of a new 
neutrino scalar force, coupled with near-gravitational 
strength to matter, can allow a resolution of the BBN 
constraints on sterile neutrinos by changing the mass pa- 
rameters in the early universe so as to prevent oscilla- 
tions. Wc shall also see that the evolution of this setup 
naturally yields the late time universe as populated with 
mass eigenstate neutrinos, providing suitable initial con- 
ditions for neutrino dark energy. 



II. NEUTRINO SCALAR FORCES 

The simplest means of generating a new force for neu- 
trinos is by the inclusion of a new, singlet neutrino and 
an associated Yukawa force. 



C D mr>vn + XAnn + Vq(A), 



(1) 



where v is the standard model left-handed neutrino, n is 
a sterile neutrino, and A is the scalar acceleron field. 

If the relic neutrinos are in the light mass eigenstate, 
then there is a contribution to the effective potential of 
the acceleron at late times 
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for a thermal background, both of which will drive A 
to large values, in order to minimize the energy of the 
neutrino background. 

However, in the early universe, thermal scatterings 
produce interaction eigenstates, not mass eigenstates, as 
long as the sterile state is light enough to be produced 
(i.e., XA = m„ < T). We must ask the question: as- 
suming the background neutrinos are presently in mass 
eigenstates, at what temperature would the sterile state 
become thermally accessible? This was studied in ^(j 
in the case that V {A) = ^ i \og{A/ Aq) and there it was 
found that XA(T) > T. However, we would like to ex- 
plore theories with quadratic potentials, both because 
quadratic potentials are very generic, and because this is 
the form of the potential in the recently proposed hybrid 
models of neutrino dark energy. 

Thus, we take a potential Vb(-4) = m 2 A A 2 /2, and find 
the background value of the acceleron field 



XA 



XT rn D 
12i74' 



777.4 



(4) 



For XA > T, we find Xm 2 D /VUmA> T. Since 
to _4 ~ Xmr> by naturalness arguments |16(, the effects 
of relativistic neutrinos are generally insufficient to keep 
the sterile neutrino out of the effective theory at temper- 
atures relevant for BBN. 

In order to have a background of mass eigenstate neu- 
trinos, and to prevent the thermalization of the sterile 
states, we must have either decohering scatterings out of 
equilibrium, such as at T < MeV, or XA > T, so that the 
heavy state cannot be produced by such scatterings. As 
we have found that the latter condition cannot be main- 
tained by neutrinos alone, we are prompted to consider 
other contributions to the effective potential in the early 
universe. 



A. Matter effects and the early universe 

Neutrinos are unique, in that they can feel the effects 
of A through mixing to light states. In contrast, no other 
fcrmion in the standard model can mix with a gauge sin- 
glet, even after electroweak symmetry breaking. How- 
ever, higher dimensional couplings between A and these 
other fermions are not only allowed, they are naturally to 
be expected from Planck-scale physics. At temperatures 
near T ~ 1 MeV, only electrons are still in equilibrium, 
and so we consider couplings of the form 
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(5) 



Such a term would arise from the higher dimension 
operator (3 e y e Aheie r / 'Mp\ where h is the standard model 
Higgs field. This is precisely the sort of operator consid- 
ered in [Til. ll8L IT^ I , which lead to new matter effects in 
neutrino oscillations. 



One might think that interactions at this strength 
(comparable to gravitational), may not have significant 
effects on the phenomenology of neutrino oscillations in 
the early universe. However, due to the enhancement of 
the high temperature, there can easily be an important 
correction to the free energy of the acceleron field, 



m e (A) 
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where 



I+{z) = dyy 2 log(l + e 

Jo 



-y/v 2 +z 



(6) 



(7) 



and g e — 4 is the number of degrees of freedom in a Dirac 
electron. At high temperatures (T 3> 77i e ), this simplifies 
considerably, yielding a total potential 



The minimum of the thermal effective potential is at 
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(9) 

To satisfy BBN constraints, we need to ensure that 
only one state is thermalized. As we have previously 
stated, this is most straightforwardly achieved by the re- 
quirement that the heavy state is thermally inaccessible 
when electron annihilation to neutrinos is still in equilib- 
rium, i.e., XA(T ann ) > T ann . If we additionally want to 
produce the proper initial conditions for neutrino dark 
energy, the heavy state must be inaccessible until after 
all scatterings have ceased, that is, XA(T scat ) > T scat , 
which results in the requirement that 
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The annihilation process v a Va e + e freezes out at 
T e w 3.2MeV and T AlT w 5.3 22}. The scattering process 
v a e «-> v a e freezes out at 1.4MeV. 

However, the coupling to electrons cannot be arbitrar- 
ily large. Such a field generates a Yukawa potential be- 
tween two atoms of the form 



V(r) 



(11) 



(3 e must be small enough to be consistent with tests of the 
gravitational inverse square law (ISL) and equivalence 
principle violation limits, as the force couples to lepton 
number L (or, equivalently, Z). These limits are typically 
quoted in terms of the parameter a, which in our system 
is 



0l< 
Anm 2 



(12) 
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where m u is the atomic mass unit. 

For a force in the range m_4 ~ 10 _7 eV, equivalence 
principle tests yield the tightest constraint, a < 10~ 6 5 
or 



/3e<4 



(13) 



This limit is roughly constant with masses down to about 
10~ 12 eV. Slightly heavier accelerons (m_4 > 10 _6 eV) 
have weaker limits on /3 e , which are ultimately domi- 
nated by ISL experiments at above ~ 10~ 5 eV. Since 
the A expectation value in the early universe drops like 
m 2 ^, these heavier accelerons tend to have significantly 
smaller expectation values in spite of the larger allowed 
values of (3 e . Note that we have ignored the possibility 
of chameleon effects [2^, 13 HH i which could weaken the 
limits. Regardless, within these constraints, eqn. 1101 is 
easily satisfied for T ~ MeV. 

Hence we see that even very weak (Planck-scale) in- 
teractions with electrons can drive the acceleron to large 
values, preventing the thermalization of sterile neutrinos 
in the early universe, and leaving the background neutri- 
nos as mass eigenstates at late times. 



B. Naturalness 

The recent SUSY hybrid models have given a set of 
models where the neutrino dark energy has all mass pa- 
rameters of their technically natural size ^(|. It is im- 
portant to consider whether the scenario described here 
can be satisfied in the context of a natural model. 

In these models, one can calculate the quantum cor- 
rections to the acceleron mass explicitly. Because the 
loops involve light fields, the masses run all the way 
from the weak scale to the neutrino mass scale, where 



the natural size of the acceleron mass is 



A 2 



At these intermediate temperatures, the natural value is 
A TTOjj/3. Expressing the previous limit with the 



nr 



replacement m A — A m D /3, the earlier limit reads 



are smaller for the same value of A, and may begin to 
thermalize earlier. 

Considering for instance the model of reference (lt| . 
we have A< 10~ 55 , m D ~ 10~ 15 eV and m A ~ 10~ 7 cV. 
The sterile state can come into thermal equilibrium at 
T ~ 70MeV for j3 e = 4. 

Once the sterile state can be produced, we must use 
a density matrix formalism to describe the dynamics, 
where 



Paa Pat 
Psa Pst 



(15) 



If p ss begins to be populated, it contributes to the ef- 
fective potential of the acceleron an amount p ss X 2 T 2 /24. 
Hence, even for small values of p ss , this can dominate 
over rn 2 ^ in the potential, driving A to smaller values, ul- 
timately allowing all sterile states to become populated. 
Therefore, we should ensure that this cannot happen un- 
til after T — 5.34 MeV, when /i, r-neutrino production 
freezes out. 

There are a number of simple possibilities to achieve 
this. First, if m A is tuned at the few percent level, 
we could accommodate A ~ 10~ 4 for run ~ 10 -15 , 
which addresses the problem. A second possibility is 
just that both A and m A are smaller than what satu- 
rates the model in [lg. For instance, one could take 



m D 



10" 2 eV, A 



10- 



10 eV (evaluated at 



the neutrino mass scale) and (3 e = 1. 

Lastly, we note one additional interesting possibility. 
Although muons become non-relativistic much earlier 
than BBN, their couplings to the acceleron are far less 
constrained, and can have a significant effect even at tem- 
peratures as low as T ~ 5 MeV. If the muons have a 
non-renormalizable coupling ATrhfj,i[j, r A/M#, then 



dV _ 2T i dI + (m^(A)/T) 
dA ir 2 dA 
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Notice that this expression can be satisfied for any m d 
by taking A sufficiently small, ~ 10~ 4 eV was shown 
to be of the appropriate size to explain the current dark 
energy. 



1. Multiple Sterile Neutrinos 

One can incorporate the masses of the heavier neutri- 
nos into these theories by the inclusion of additional ster- 
ile states with somewhat larger Dirac masses. However, 
continuing with the naturalness arguments, these neutri- 
nos with larger Dirac masses should have smaller cou- 
plings to the acceleron. As a consequence, their masses 



which yields 

A(T = 5.34 MeV) = 10* MeV (^) " (H^X) . 

^ (17) 

This would prevent the sterile state mixing with [i and r 
neutrinos from being populated until after v^ T produc- 
tion had ceased. The couplings to the electron, however, 
would still be essential to prevent the sterile states mix- 
ing strongly with v e from thermalizing. Such a model 
could arise in extra dimensional scenarios, for instance, 
where the acceleron was dominantly localized on the 
muon brane. 

Ultimately, if we are to consider multiple sterile neutri- 
nos instead of just one light one, we are forced to consider 
as well certain tensions. However, the tensions are not 
severe, and straightforward solutions exist. 
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C. Initial Conditions and Acceleron Dark Matter 

We have assumed up to this point that the acceleron 
lies at the minimum of its potential, but it could have 
begun with energy in the form of coherent oscillations. 
There are a number issues we must consider. The first 
is that of the usual moduli problem, that the energy in 
acceleron oscillations not cause early matter domination. 
This implies 

^<A 2 )(^f) 3 </w, (18) 

where Ti, the temperature where A begins to oscillate, 
is defined by the condition m A = iJ(Ti), and A% is the 
initial acceleron vev, which we assume can be very large. 
Eqn. El then yields the condition 

(A><(10 12 GeV) 2 fH^£XV /2 , (19) 

V m A J 

so that for a wide range of Ai, the energy in acceleron 
oscillations does not exceed the total energy in dark mat- 
ter. 

A second concern is that prior to BBN, when the neu- 
trino scattering and annihilation processes occur at a 
rate 7 which exceeds the oscillation rate of the acceleron 
about its minimum, 7 > m A , it is possible that for cer- 
tain periods of time we will have m n < T, allowing the 
sterile states to be thermalized. Prior to electroweak 
symmetry breaking (EWSB), the Yukawa coupling be- 
tween the sterile neutrino, active neutrino and Higgs is 
too small to allow thermalization. Hence we are only wor- 
ried about thermalization during oscillation after EWSB. 
Should this occur, there can be large contributions to the 
acceleron effective potential, driving (A) towards smaller 
values, allowing complete thermalization to occur. 

The amount of time that the oscillating field stays be- 
low XA < T for (.A 2 ) (A) 2 can be estimated as 

T T 

XA ~ Xm a (A 2 ) 1/2 (20) 

The scattering rate is ~ lCr 22 (T/MeV) 5 . Hence, for 
an initial vev of 10 12 GeV or less, we see that it is in- 
evitable that there will be many scatters during the pe- 
riod when the sterile state can be thermalized. 

We have already seen that the finite temperature ex- 
pectation value for A falls like T 2 , while the amplitude 
of oscillations (A 2 ) falls like T 3 , so this requirement be- 
comes more stringent at later times. Once m A > t _1 we 
can just use the average value of (A) , since there is insuf- 
ficient time to thermalize. Requiring that the amplitude 
of oscillation is smaller than the expectation value at T x 
when r _1 = m A gives us an initial amplitude 

x 1 1 mm\M 2 x e V m A J 

(21) 



Satisfying this constraint ensures that there is no con- 
flict from BBN. However, if we additionally want to en- 
sure that the oscillations do not spoil the initial condi- 
tions for neutrino dark energy, we must concern ourselves 
with what occurs at somewhat later times. After BBN, 
electrons become non-relativistic and quickly lose any sig- 
nificant contribution to the effective potential. At this 
point, only the neutrino contributions to the acceleron 
potential are relevant. 

If the acceleron can oscillate all the way to zero from 
this minimum, then it is possible that the neutrinos hop 
from one mass eigenstate to the other when they are de- 
generate (i.e., when A = 0). This is particularly possible 
because the mixing terms go to zero more slowly than 
the difference of the diagonal eigenvalues, leading to a 
particularly non-adiabatic situation. 

At small values for A, the contribution to the potential 
from the neutrinos reaches a maximum of T 2 m 2 D /2A. If 
there is less energy in the oscillations at this point than 
this, the acceleron will never oscillate to zero. Thus, we 
have a limit on the initial amplitude 

^>".i<A 2 >(|) S . (22) 
Taking T e = m e , we find 

(«<.C)"(«v)'«' e» 

III. NON-STANDARD MATTER EFFECTS 

The lower limit on /3 e indicated in ean. 1101 gives rise to 
the possibility of experimentally verifiable signals in neu- 
trino oscillation experiments. In it was shown that 
these new matter effects, for typical earth densities could 
overwhelm the traditional MSW effect in earth and give 
rise to significant signals in neutrino oscillation experi- 
ments. The new matter effects give rise to an effective 
potential 

V{A) = ^-p & A+\m\A\ (24) 

where p e ~ 10~ 3 g/cm 3 is the typical electron density in 
earth. 

Then the sterile neutrino mass will change in earth as 

/ 10 _6 eV\ 2 

Am„ = AA.4 ~ 10 _1 eV xAft . (25) 

V m A J 

The requirement that mass eigenstates be populated to 
produce dark energy at late times, however, gives a lower 
limit on the product \/3 e /m A through ean. 1101 Combin- 
ing with eqn. ^| (recognizing that we must evaluate the 
mass here at zero temperature), we find a minimal ex- 
pected size of matter effects which is 

Am„ > 10 _1 eV. (26) 
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If the acceleron also couples to baryons, the effect could 
be much stronger. Since forces couping to baryon number 
do not violate the equivalence principle as strongly as 
those coupling to electron number, the effects from as 
could be roughly fifty times larger. We emphasize that 
the above value is merely a lower limit on this effect. 

Thus one expects matter effects to have significant im- 
pact in experiments which may constrain the sterile neu- 
trino mass and mixings, for example LSND [2(|, KAR- 
MEN H3 and MiniBooNE [2^. It has been show that 
such sensitive matter effects can help to reconcile the 
LSND result with the Bugey 0, CDHS HfJ and KAR- 
MEN experiments, where the combination of limits from 
the experiments disfavors 3+1 and 2+2 fits to the neu- 
trino oscillation data [19| . In addition, such large matter 
effects may modify MiniBooNE 's sensitivity when com- 
pared with LSND. 

We expect smaller, though perhaps not insignificant, 
effects on the mostly active mass eigenstates: 



Am„ = Ar 



m 



D 



n m n (m n + Am n ) ' 



(27) 



where m n is the vacuum sterile neutrino mass. 



A. Consequences for future neutrino experiments 

We have shown that requiring a set of initial conditions 
which satisfy BBN constraints and give the appropriate 
initial conditions for dark energy at late time give rise 
naturally to non-standard matter effects of a size which 
could be detectable at upcoming experiments. 

One interesting possibility, noted elsewhere [Isl H^ . 
is that such a matter effect might serve to explain the 
discrepancy between LSND and the null short-baseline 
experiments. In particular, because the pathlcngth at the 
Bugey experiment is mostly air, while LSND is mostly 
earth, the allowed parameter space for the LSND signal 
is significantly widened. This would then allow a signal 
at MiniBooNE, even in a region of the parameter space 
which is already excluded in the standard 3+1 scenario. 
Such a signal might reasonably not agree with the LSND 
results, if the intervening material is sufficiently different. 

Future short baseline experiments also hold the possi- 
bility to detect this signal. In short baseline experiments, 
the amplitude for oscillation goes as 

1 - P ee = sin 2 20 sin 2 {Am 2 L/AE) ~ ( 28 ) 



.016 xf^) V ™" 
KleV) VIO" 1 



eV 



50 m 



4E 2 
)0 Me 



L \ 2 /l00MeV\ 2 



where the second equality holds in the limit that mjj <C 
m n and L < AE / m 2 . Hence we see that should m n go up 
in matter, the effect would be to increase the oscillation 



probability into sterile states. An experiment with the 
sensitivity of MiniBooNE may detect such an effect with 
mj ~ 0.2 eV and Sm n ~ 1 eV. Alternatively, a two- 
detector experiment to search for #13 , might perform two 
high statistics runs, with different amounts of material 
between the neutrino source and the near detector, and 
thus place limits on this effect. 



IV. SUMMARY 

The past decade has been extremely exciting for be- 
yond the standard model physics, in that we have now 
found convincing evidence for neutrino mass and dark 
energy. A possible connection between the two requires 
the introduction of a singlet neutrino and a new scalar 
field, and thus forces us to consider their properties in 
the early universe. 

Couplings of the scalar to electrons with gravitational 
strength are to be naturally expected. We have shown 
such interactions significantly alter the early universe cos- 
mology of the system. In particular, the potential cre- 
ated by this coupling generally forces the singlet neutrino 
to be sufficiently heavy that it is thermally inaccessible 
through BBN, when neutrino scattering and annihilation 
processes which can populate the sterile state decouple. 

A sufficiently large coupling would naturally lead to 
significant effects at terrestrial neutrino experiments, 
particularly with regard to the mass of the singlet neu- 
trino. Such effects change the relative sensitivity of dif- 
ferent experiments depending on the electron density of 
the medium of neutrino propagation. 

The process by which the singlet state is thermally in- 
accessibly at early times also leaves the relic neutrinos 
in mass eigenstates at late times. This is important as 
it provides the proper initial conditions for recently con- 
sidered "hybrid" models of neutrino dark energy. Such 
an initial condition is generic with these couplings, and 
should be a motivation to consider new matter effects a 
natural prediction of neutrino dark energy theories, pos- 
sibly detectable at MiniBooNE or a double-CHOOZ style 
experiment. 

Note added: As this work was to be submitted, [3l| ap- 
peared, which considers the consequences of MiniBooNE 
for MaVaN scenarios. 
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